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Abstract 

We present numerical evidence that there are two distinct universality classes 
characterizing driven interface roughening in the presence of quenched disor- 
der. The evidence is based on the behavior of A, the coefficient of the nonlin- 
ear term in the growth equation. Specifically, for three of the models studied, 

A — ► oo at the depinning transition, while for the two other models, A — > 0. 
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The motion of a nonequilibrium interface in a disordered environment has attracted 
much attention. Fluid flow in a porous medium is a typical experimental realization of these 
phenomena, but applications range from wetting phenomena to the motion of flux lines in 
the presence of disorder The origin of different universality classes is well understood 
for growth in which thermal (or time-dependent) noise dominates the roughening process. 
However, recently several independent studies Jl| ^] have noted that quenched noise — which 
is independent of time and depends only on the position of the interface — may change the 
universality class, generating interfaces with anomalously large roughness exponents. 

In the typical d- dimensional interface characterized by a height h(x,t) moves in 

a (d + l)-dimensional disordered medium. The randomness of the medium can be described 
by a quenched noise r](x,h). In the presence of an external driving force F, the simplest 
growth equation describing the zero-temperature dynamics of the interface is [fJJ 

d t h = F + uV 2 h + r](x,h). (1) 

The z/V 2 /i term mimics a surface tension and acts to smooth the interface, while the quenched 
noise r)(x, h) works to roughen the interface. It is generally assumed that the quenched noise 
has zero mean and is uncorrelated. 

An interface characterized by ([!]) moves with a finite velocity vq if the driving force 
exceeds a critical value F c , while for F < F c it is pinned by the disorder. When F —>■ F c , 
one finds 

~ f, (2) 

where / = (F — F c )/F c is the reduced force and 9 is the velocity exponent. 

Recently, a number of analytical and numerical [Q-|7J] studies focused on understanding 
the nature of the depinning transition and obtaining accurate estimates for the critical 
exponents. Renormalization group (RG) studies of Eq. (|1|) find a roughness exponent 
a = (4 — d)/3, but a number of numerical models [Q-^] revealed exponents whose values 
can be quite different from the RG predictions. Wetting fluid invasion models gave a ~ 
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0.8 for (1 + 1) dimension, and investigation of the random field Ising model (RFIM) in 
(2 + 1) dimension gave a ~ 0.67 ||. Solid-on-solid type models gave a ~ 0.63 for (1 + 1) 
dimension, and a ~ 0.48 for (2 + 1) dimension @||, while a discretized solid-on-solid model 
of Eq. ([I]) gave a ~ 1.25 for (1 + 1) dimension, and a ~ 0.75 for (2 + 1) dimension J7|. A 
similar scattering is found for the values of the other exponents characterizing the depinning 
transition (F = F c ). 

Here we report simulations of five distinct models that have been introduced to inves- 
tigate the motion of an interface in the presence of quenched disorder ^-0]. Our findings 
suggest the existence of two different universality classes. One universality class is described 
by the nonlinear growth equation || 

d t h = F + uV 2 h + X(Vh) 2 +rj(x,h), (3) 

and we find that for models in this universality class 

A ~ /-*. (4) 

The second universality class is described, at the depinning transition, by (|). We propose 
that the existence of the two universality classes is the origin of the systematic differences 
between the exponents predicted by RG calculations and estimates from numerical simula- 
tions. We find that measuring A can give information about the universality classes to which 
a given growth process belongs. 

To calculate A, we follow Ref. and impose a "tilt" of slope m on the interface. For 
a (1 + l)-dimensional system, we consider a lattice with L columns, and "build in" the tilt 
by implementing helicoidal boundary conditions, h(0,t) = h(L,t) — Lm and h(L + l,t) = 
h(l,t) + Lm. For a (2 + l)-dimensional system, we use the same boundary conditions, so 
the tilt occurs only in one direction. 

For a model described by Eq. (|l|), the local velocity v is independent of the tilt. However, 



if a nonlinear term is present in addition to the linear term, then from fl3|) follows that flO 



v = vq + Am . (5) 



Thus by varying the tilt m, we can test for the presence of nonlinear terms in the growth 
equation and calculate the coefficient A. 

First, we treat the model introduced in Ref. Q, for which it was shown, in (1 + 1) 
dimension, that the interface at F c is pinned by a directed percolation (DP) cluster [fH|, 
and that the critical dynamics are controlled by a divergent correlation length parallel to the 
interface £y ~ with v\\ ~ 1.73. This model, referred to as "DP-1", excludes overhangs, 
and gives rise to a self-affine interface at the depinning transition, with a roughness exponent 
a ~ 0.63 |T[|. 

In the DP-1 model, in (1 + 1) dimension, we start from a horizontal interface at the 
bottom edge of a lattice of size L. At every site of the lattice we define a random, uncorrelated 
quenched variable, the noise T], with magnitude in the range [0, 1]. During the time evolution 
of the interface, we choose one of the L columns at random. If the difference in height 
to the lowest neighbor is larger than (+1), this lowest neighboring column grows by one 
unit. Otherwise, the chosen column grows one unit provided the noise on the site above 
the interface is smaller than the driving force F. The unit time is defined to be L growth 
attempts. 

We measure the velocity of the interface for different reduced forces / and different tilts 
m. The results for (1 + 1) dimension are shown in Fig. |l]a. For a fixed force /, we find that 
the interface velocity changes with the tilt m, indicating the existence of nonlinear terms. 
Near the depinning transition (/ — > 0), the velocity curves become "steeper" and from (|5|), 
we infer that A must increase. 

To measure A, we first attempt to fit a parabola to the tilt-dependent velocities in the 
vicinity of zero tilt. The calculations indicate that as we approach the depinning transition, 
A diverges according to Eq. (Q). However, in the vicinity of F c , the velocity curves lose their 
parabolic shape for large tilts (see Figs |l]a and|2|), indicating the presence of other terms not 
included in (pp. 

We can understand the breakdown of (^) for large m using scaling arguments. Substi- 
tuting Eqs. (|]) and ([|) into (§]), we find 
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v{m,f) oc f + af-^m 2 . 



(6) 



Equation (fj|) indicates that the velocity curves corresponding to two different forces fi 
and $2-1 with fi > f 2 , will intersect at a tilt m x (see Fig. |2|). For tilts greater than m x , 
v(m, fi) < v(m,f 2 ), a clearly unphysical prediction; the average velocity, for the same tilt, 
should be larger for the larger force. Thus the velocity cannot follow a parabola for arbitrarily 
large m, and a crossover to a different behavior than that of Eq. ([]) must occur for values 
of the tilt larger than m x . 

Letting (fi — j'2) — > 0, we find from (|6|) that the crossing point of the two corresponding 
parabolas scales as 



ml ~ f 



(?) 



Equations (J|) and (0) motivate the scaling form for the velocities 



(8) 



where g(x) ~ const, for x <C 1, and ^(x) ~ x e/(e+<t>) f or ^ 3> 1 [12|. Figure 0a shows 
the data collapse we obtain using fl8|), and the results of Fig. |TJa rescaled with exponents 
9 = 0.64 ± 0.08, = 0.64 ± 0.08 for (1 + 1) dimension 

The scaling behavior (|8]) is not limited to the DP-1 model in (1 + 1) dimension, for (2 + 1) 
dimension and for the models introduced in Refs. |||6|] we find a very similar behavior. We 
refer to these models as "DP-2" || and "Parisi" 0. We simulated them in (1 + 1) dimension, 
and were able to rescale the velocities according to (S) using the exponents presented in Table 



Another model we studied was the RFIM, which allows for overhangs; and for certain 
values of its parameters can be mapped to percolation 0. In the RFIM, spins on a square 
lattice interact through the Hamiltonian 



H=- y £S i S j - y £[F + ri(i,h)]Si 



(9) 
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where Si = ±1, F now denotes the external magnetic field, and 77 is the time- independent 
local random field (i.e., quenched noise) whose values are uniformly distributed in the interval 
[—A, A]. The strength of the quenched disorder is characterized by the parameter A. At time 
zero, all spins are "down" — except those in the first row, which are initially up. During the 
time evolution of the system, we flip any down spin that is "unstable," i.e., whenever the flip 
will lower the total energy of the system. The control parameter of the depinning transition 



is the external magnetic field F; the unit time corresponds to flipping all unstable spins |L4 

For dimension (1 + 1), there are two morphologically-different regimes, depending on 
the strength A of the disorder (i.e., of the random fields). For A > 1.0, the interface is 
self-similar (SS), while for A < 1.0 it is faceted (FA). For dimension (2 + 1), there is again a 
FA regime (A < 2.4), a SS regime (A > 3.4), and also a self-affine (SA) regime in between 
(2.4 < A < 3.4) [[3]. The SA regime, which exists only for (2 + 1) dimension, is the only 
regime of the RFIM for which either Eqs. ([I]) or @ could apply. In the SS regime, the 
interface is not single-valued, while in the FA regime, lattice effects dominate the growth. 

Our results show that for the FA regime, the RFIM behaves in a similar fashion to the 
other three models, in that the coefficient of the nonlinear term diverges at the depinning 
transition. However, although (g) is still valid for the SA and SS regimes, we find a negative 
4>, thus A — > 0. This behavior can be understood, for the SS regime, by considering that 
near the depinning transition, the morphology of the interface corresponds to the hull of a 
percolation cluster, which has no well-defined orientation 0. Thus a change in the boundary 
conditions will not affect the growth process, and we cannot expect any divergence of a 
possible nonlinear term when the magnetic field approaches its critical value. On the other 
hand, for large fields, the effect of the quenched disorder diminishes, and we can observe an 
average interface orientation. For such values of field, we expect the presence of nonlinear 
terms to be felt. Although for the SA regime the behavior of A is similar (see Figs. |I]b and 
^|b), the reasons so far cannot be understood. 

These results lead us to conclude that in the SA regime the RFIM belongs to the univer- 
sality class of Eq. ([!]). This conclusion is further supported by the agreement between the 



numerically determined exponents, a ~ 0.67 and 9 ~ 0.60 for (2 + 1) dimension, and the 
RG predictions, a = 9 = 2/3 §. 

Finally we studied the discretized solid-on-solid version of Eq. ([!]), referred to as "SOS-1" 
0], and find that for any reduced force, A = 0. 

The results of Table | show, for (1 + 1) dimension, a clear separation into two groups in 
the values of the critical exponents for the five models studied. In the following we argue 
that this separation reflects the existence of two distinct universality classes, described by 
the two continuum growth equations, ([!]) and (§). For the SOS-1 model and for the RFIM, 
in the SA regime, we find that A either is zero or goes to zero at the depinning transition. 
Thus the scaling behavior of these models should be correctly described by (|]). For the 
DP-1, DP-2 and Parisi models we observe a divergent A, indicating that nonlinearities are 
relevant near the depinning transition. Thus to properly describe the scaling properties of 
these models it is necessary to study (||), since ([[]) does not includes the nonlinear term 
X(Vh) 2 . Further evidence of the existence of the two universality classes is given by the 
values of roughness exponents. The models for which A diverges at the depinning transition 
predict a ~ 0.63, in agreement with the mapping to DP. On the other hand, models 
in the universality class of Eq. ([[]) 0,0], gave roughness exponents typically larger, in better 
agreement with the RG predictions ||. Finally, we propose the study of the behavior of A 
at the depinning transition as a general method for identifying the universality class of a 
given growth process in disordered media. 

This work is part of the Ph.D. thesis of L. A. N. Amaral, to be submitted to Boston 
University. We thank R. Cuerno, S. Havlin, H. Makse, O. Narayan, T. Nattermann and 
especially S. V. Buldyrev for useful suggestions. L. A. N. Amaral acknowledges a scholar- 
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FIGURES 

FIG. 1. Dependence on the tilt m of the average velocity, (a) in the DP-1 model and (b) in the 
RFIM. Data for different forces / are indicated by different symbols, ranging from 0.016 (bottom 
curve) to 0.350 (top curve) for the DP-1 model, and from 0.014 (bottom curve) to 0.143 (top 
curve) for the RFIM. In (a) we show velocities, for the DP-1 model, for a system of size 512 in 
(1 + 1) dimension. In (b) are plotted the velocities for the RFIM in the SA regime (A = 3), for a 
(2 + l)-dimensional system of size 40 x 40. 

FIG. 2. Here we exemplify the "noncrossing" effect on the velocity parabolas. We show a perfect 
parabolic behavior for two different forces, /i > ji (dashed lines) as predicted by Eq. (5). Also 
shown is the "curving back" of the velocity curve for the smaller force fa (solid line) in order not 
to cross the velocity curve for f\. 

FIG. 3. Data collapse according to (8), using the same symbols for the velocities shown in Fig. 
1. In (a) we present the rescaled results for the DP-1 model in (1 + 1) dimension and in (b) the 
rescaling of the (2 + l)-dimensional results for the RFIM, in the SA regime (A = 3.0). 
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TABLES 

TABLE I. Exponents for the five studied models (see definitions in the text). A negative value 
of <j> means that A — > when / — ► 0. We argue in the text that the models above the horizontal 
line (DP-1, DP-2, and Parisi) belong to the universality class of Eq. (3) and can be mapped, in 
(1 + 1) dimension, to DP. The models below the line belong to the universality class of Eq. (1). 



Model 




fl + 1) 

9 


dimension 


(2 + 1) 

9 


dimension 


DP-1 




0.64 ±0.08 


0.64 ±0.08 


0.80 ±0.12 


0.30 ±0.12 


DP-2 




0.59 ±0.12 


0.55 ±0.12 






Parisi 




0.70 ±0.12 


0.65 ±0.12 






RFIM 


SA 
SS 


0.31 ±0.08 


-0.65 ±0.13 


0.60 ±0.11 


-0.70 ±0.11 


SOS-1 




0.26 ±0.07 
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